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The problem of the metastable phase decay remains one of the most 
actual questions in the field of kinetics of the first order phase transitions. 
The first theories 0, 0] determined only the characteristic time lag of the 
nucleation process. The first well grounded theory which allows to describe 
the evolution of the system during the nucleation period and to find the 
form of the droplets size spectrum was presented in [Q] where the iteration 
procedure was formulated. 

Now the theoretical description is going to be spread on the systems 
with the several types of heterogeneous centers. Unfortunately the direct 
generalization of the iteration method failed due to existence of the specific 
interaction between different sorts of heterogeneous centers through the vapor 
consumption. 

The first attempt to solve this problem was done in where the structure 
of announced description was composed by many asymptotic cases including 
the intermediate case where the new special monodisperse approximation 
was suggested. 

This monodisperse approximation was well grounded but there was no 
evident correspondence between the iteration procedure and the results pre- 
sented in [0. Following the suggestion of F.M.Kuni to seek for this corre- 
spondence we shall present here the iteration solution of this problem. 

We shall consider the physical situation from and use all definitions 
introduced there. 

The structure of the publication will be the following: 
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• In the first part the iteration procedure will be presented. 

• In the second part we shall calculate iterations on the base of the prop- 
erty of the "avalanche consumption" of the metastable phase. 

• In the third part we shall calculate the iterations based on the special 
monodisperse approximation. 

1 Iteration method 

We shall present all formulas for the case of two sorts of heterogeneous cen- 
ters. They will be marked by two subscripts A and B. The generalization 
for the arbitrary number of heterogeneous centers types is evident. It is one 
of the real advantages of the presented theory. 

The system of the condensation equations according to 0] can be pre- 
sented as following 

Ga{z) = Fa j\z - xfeM^Aia^) - ^))eA{x)dx 
Jo 

Gsiz) = Fb I\z - xf exp(rB(C(x) - ^))eB{x)dx 
Jo 

^ = C{z) + Ga{z)+Gb{z) 

eA{z) = exp{-KA rexp{TA{ax) - <!>))dx) 
Jo 

eB{z) = exp{-KB I exp(rB(C(x) - ^))dx) 
Jo 

with six parameters Fa, Fb, Ka, Kb, Ta, F^. Note that the definition of 
Ta, slightly differs from the standard one - it is $ times less. 

The presented system can be directly seen from the condensation equa- 
tions system for one type of heterogeneous centers and there is no need to 
discuss it more. The problem is how to solve it. 

According to the scale invariance one can put Fa = 1, Ta = I- According 
to the symmetry of the types one can choose them as to have Fb < 1. These 
simplifications will be interesting only in numerical modeling and here we 
shall keep all parameters. 
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We shall define the iteration procedure by the following relations 

Gm+i{z) = Fa I\z - x)3exp(r^(0(a;) - ^))eM{x)dx 
Jo 

GBt+i{z) =Fb f {z- xfexp{rB{Ci{x) - <^))eBi{x)dx 
Jo 

$ = Ci{z) + GAiiz) + GBiiz) 

OAi+i{z) = exp{~KA f exp{TA{Ci{x) - ^))dx) 
Jo 

OBi+i{z) = exp{-KB f exp{TB{Ci{x) - ^))dx) 

Jo 

The number of the iteration approximation is marked by the ordinary sub- 
script instead of the capital letter which marks below the sort of the centers. 
The values without the number of iteration will mark the real solutions. 
The initial approximations are the following 

Gao = 0, Gbo = 0, 6ao = 1; (^BO = 1 
Then one can note the following important chains 



Gao < Ga2 < ■ 


•• < GA2i < ■■ 


.<Ga<. 


■■ < GA2i+l < 


••• < Ga3 < Gai 


Gbo < Gb2 < • 


■■ < GB2i < ■■ 


. <Gb < ■ 


.. < GB2i+l < 


... < Gb3 < Gbi 


^AO > dA2 > 


■■■ > dA2i > ■■ 


..>eA>. 


■■ > 0A2i+l > ■ 


■ ■ > Qaz > Sai 


Gbo > dB2 > 


■■■ > dB2i > ■■ 


.>eB>. 


■■ > dB2i+l > ■ 


.. > 9b3 > Obi 


Co > C2 


> ••• > C2i > ■ 


..>(>.. 


■ > (2i+l > ••• 


> Cs > Ci 



for all values of arguments. 

These chains prove the convergence of iterations and allow to estimate 
the accuracy at every step of the calculation. 

The problem is to calculate the iterations. In the situation with one type 
of heterogeneous centers the system of condensation equations can be gotten 
if we chancel Gb and the index A. The same thing we have to do with the 
iteration procedure and one can get the mentioned properties. In fl]] only 
one iteration for G and two first iterations for 6 were calculated. The further 
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iterations can not be calculated analytically. In the situation of nucleation 
in the system with one type of heterogeneous centers it is sufficient. But in 
the situation with several types of heterogeneous centers it isn't so. To show 
this we has to recall the physical sense of the first iterations. 

The first iteration for G is calculated on the base of unexhausted number 
of heterogeneous centers. So, the behavior of C, is wrong when the power of 
exhaustion is essential. But when the centers are exhausted the number of 
droplets is known - it is equal to the initial number of centers. This primitive 
notation is the reason why the second iteration for 9 gives already suitable 
result. 

In the situation with the several types of heterogeneous centers the situa- 
tion is different. Really, it is quite possible to have the exhaustion of the first 
type of heterogeneous centers and the moderate exhaustion of the second 
type centers. Then the exhaustion of the first type centers has to be taken 
into account in calculation of the supersaturation which is necessary to get 
the number of the droplet appeared on the second type centers. But as far 
as the power of exhaustion of the second type centers is moderate one can 
not say that the number of the droplets on the second type centers is equal 
to the total number of the second type centers. So, we have to know the 
behavior of supersaturation in all cases and it is impossible to do already in 
the first iterations. Then we have to calculate further iterations. 

Unfortunately it is impossible to calculate the further iteration approx- 
imations without any simplifications. That's why two ways of possible ap- 
proximations are presented in the next sections. 

2 Avalanche consumption of the metastable 
phase 

The sequential calculation of iterations gives 

Gai ~ ^a-^ 
9ai = exp{-KAz) 
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6bi = exp{-KBz) 

Ga2 = Fa I\z - xf exp(-rA(F4 + Fb)x^/A) exp{-KAx)dx 
Jo 

Gb2 = Fb I\z - xf exp(-rB(FA + Fb)x^/4) exp{-KBx)dx 
Jo 

The problem is how to calculate the integrals for 6*^2! Gb2- To calculate 
these integrals we shall extract functions 

LfA = exp(-r^(G^ + Gb)) 

cpB = exp(-rB(G'^ + Gb)) 
These functions can be calculated in the iteration approximation as 

ipAi = exp{-TA{GAi + Gbi)) 

ipBi = exp{-TB{GAi + GBi)) 

One can see that ipAi, fBi can be well approximated as the step functions. 
Really, exp(— x'') can be approximately interpreted as the step function^. 
Then 

Ga2 = Fa (z — x)^ exp(—KAx)Q{ziA — x)dx 
Jo 



Gb2 = Fb / {z — x)^ exp{—KBx)Q{ziB — x)dx 
Jo 



where 
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2^171 = 



ZlB = 



TAiFA + FB) 

4 



)l/4 

a/4 



■Tb{Fa + Fb) 

and these integrals can be easily calculated]^. 
For z > Z\A we have 

3 

Ga2 = ^PAiZ' 
i=0 



"'^This was the reason why in experiments one can observe characteristic time lag and 
this characteristic was introduced in the theory of the metastable phase decay. Then the 
first theoretical descriptions adopted this fact as the given one without any justification. 

^More correctly one gas to multiply zia and zib on In"'^^'* 2. This corresponds to the 
definition of the halfwidth at the halfheight. 



5 



where 

su-iy fziA 

PAi^^TTT. ^\Pa X ex.Yi{-KAx)dx 

i\{6 — ly. Jo 

are some constants which can be calculated analytically. 
For z > ZiB we have 

3 

Gb2 = ^PBiZ' 
i=0 

where 

su-iy fziB 

PBi^TTT;: t^Fb X 'ex.p{-KBx)dx 

z!(3 — z)! Jo 

are some constants calculated analytically. 
For z < Z\A we have 

X 

Ga2 = FaQ-^v{-'^az)^{{zKj^'^q-x:p{zKa) - 
ZizK a)^ %yi\>{zK a) + ^zKAQyiv{zKA) — 6(exp(ir^2;) — 1 
For z < ziB we have 

1 



Gb2 = FBCxp{-KBz)-^[{zKBf expizKe) - 
3{zKBy exp{zKB) + 6zKb exp{zKB) — 6{exp{KBz) — 1 
The calculation of 9a2, Oa2 can be done in the same way 

eA2 = ^M-Ka r eM-^A{FA + FB)'^)dx) 
Jo 4 

Ob2 = bM-Kb r eM-^B{FA + FB)^)dx) 
Jo 4 

Then 

eA2 = eX.Y>{-KA[Q{ZiA - Z)z + Q{Z - Zia)Zia]) 

9b2 = exp{-KB[Q{ziB - z)z + Q{z - zib)zib]) 
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Now wc have to calculate 9a^^ Ob3- One can note that functions <fA2, ^B2 
can be also presented as the step functions. Namely, all functions 

exp{-rA{FAexp{-KAz)^[{zKAfexp{zKA) - 
3{zKAf exp{zKA) + 6zKAexp{zKA) - 6iexp{KAz) - 1)] + 
FBexp{-KBz)^[{zKB)^eMzKB) - SizKB^ eM^^B) + 

GzKBeMzKB) - 6{exp{KBz) - 1)])) 
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exp{-TA{FAexp{-KAz)—^[{zKAyexp{zKA) - 3{zKAy exp{zKA) + 

3 

&zKAex^{zKA) - Q{ex^{KAz) - 1)] + Y.PBiz")) 

1=0 

3 3 

exp(-r^(^p^i^' + Y^pBiz')) 

1=0 1=0 

3 1 

exY>{-VA{Y.PAiz' + FBQ-x.Y>{-KBz)—^[{zKBfQ-x.^{zKB) - 

i=0 

^izKsf expizKe) + 6zKb exp{zKB) - 6{exp{KBz) - 1)])) 
exp(-rB(FA exp{-KAz)-^[{zKAf exp{zKA) - 
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3{zKAfexp{zKA) + Q>zKAex^{zKA) - ^{exp{KAz) - 1)] + 

FBexp{-KBz)^[{zKBfexp{zKB) - SizKsf cxp{zKb) + 

6zKb exp{zKB) — 6{exp{KBz) — 1 



exp{-rBiFAexp{-KAz)^[{zKAfexp{zKA) - SizKAf exp{zKA) + 

3 

6zKAexp{zKA) - 6(exp(X^z) - 1)] +J2PBiz')) 

1=0 



K 



A 



3 3 

1=0 1=0 



e-x.^{-TB(^PAiz' + FBe-x.^{-KBz)—^[{zKBfe-x.v{zKB) - 

i=0 B 

SizKB^eM^^B) + 6zKBexp(zKB) - 6(exp(XB^) - 1)])) 

have the step-hke behavior. This is the central point of our calculations. 
This property can be directly seen by calculations. 

The possibihty to see the step-like behavior of exponent of the supersat- 
uration deviation directly on the base of explicit expressions is the evident 
advantage of the iteration method. Otherwise we have to prove this prop- 
erty in the general situation taking into account that exponent of the super- 
saturation deviation lies between exp{—z^) and exp{—z^) after the suitable 
r enormalization . 

Now to calculate 9a3, Obs we have to define the values Z2a, Z2b by equa- 
tions 

^a{Ga2{z2a) + GB2iz2A)) = 1 
^b{Ga2{z2b) +Gb2{z2b)) = 1 

Then 

61^3 = exp{-KA[Q{z2A - Z)Z + @{Z - Z2a)z2a]) 

9b3 = ex.p{-KB[Q{z2B - z)z Q{z - Z2b)z2b\) 
The calculation of and G^b can be done by the following way 

Ga3 = Fa f\z - x)^exp{-TA{GA2 + Gb2)) 
Jo 

exp{-KA[Q{ziA — x)x + Q{x - ziA)ziA])dx 

Gas^Fa (\z-xfQ{z2A-x) 
Jo 

exY){—KA[Q{ziA — x)x Q{x — ziA)ziA\)dx 

Gb3 = Fb I\z - xf exp{-TB{GA2 + ^^2)) 
Jo 

exp{-KB[Q{ziB - x)x + Q{x - ZiB)ziB\)dx 
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Gbz = Fb j\z - xfQ{z2B - x) 
Jo 

exp{—KB[0{ziB — x)x + Q{x — ziB)ziB])dx 
These integrals can be easily taken in analytical form. For Gas one can 

get 

• For z > Z2A with arbitrary zia we have 

3 

Ga3 = Yl ^'PAi 
1=0 



where 

PAi = ... Fa / x'[exp{-KAZiA)0{x - zia) + 

iU6 — i)l Jo 



x'[e5 

exp{—KAx)Q{ziA — x)]dx 



can be calculated analytically. 

• For z < ■min{ziA, Z2a} 

Gas ^ Fa [ {z - xf ex.p{-KAx)dx 
Jo 

has been already calculated. 

• For ZiA < Z < Z2A 

Gas = ^'PAi + Fa — exp(-i^A^iA) 

i=o 

where 



Pai^Fa—- ^/ x'exp{-KAx)dx 

i\{6 — i)\ Jo 



All PAi are constants and can be calculated analytically. 
For Gbs one can get 
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For z > Z2B with arbitrary Zib 



3 

Gb3 = z'pBi 
i=0 



where 



PBi = -777. t-.Fb \ x'[exp(-irB2;iB)©(x - Z^b) + 

— l)\ Jo 

exp{—KBx)Q{ziB — x)]dx 

• For z < ■min{ziB, Z2b} 

Gb3 Fb [ (z- eKp(-KBx)dx 
Jo 

has been already calculated. 

• For ZiB < Z < Z2B 

Gb3 = X! ^'PBi + Pb exp{-KBZiB) 



i=0 

where 



PBi = Fb7777 ^ / x' fi:^-^{-KBx)dx 

l\(6 — l]\ Jo 



■i!(3-i)! Jo 

All PBi are constants and can be calculated analytically. 

One can easily see that zia < Z2a, zib < Z2b and the third opportunity 
can take place. 

Moreover, if we define ZiA, zib by equalities 

^A{GAi{ZiA) + GBi{ZiA)) — 1 

^ b{G Ai{ZiB) + GBiiZis)) = 1 

then we can get we following chains of inequalities 

ZlA < ZsA < ■■■ < Z2i+IA < ■■■ < Za< ■■■ < ^2iA < ■■■ < ^4A < 
ZlB < Z3B < ■■■ < Z2i+IB < ■■■ < Zb < ■■■ < Z2iB < ■■■ < Z4B < Z2B 
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where za, zb are defined by 

rA(GA(zA)+GB(zA))^l 
rB{GA{zB)+GB{zB)) = l 

Then we can calculate Ga4: Gba- In the initial expressions 

Ga4. = Fa j\z - xf&c^{-VA{GAz + GBz)eAzdx 
Jo 

Gba = Fb [\z - x)^ exp(-r^(G'A3 + GBz)eB^dx 
Jo 

we have to use two characteristic lengths Z2 and z^. Then 

Gaa = Fa f e{zsA -x){z- x)^ 
Jo 

[exp{-KAx)Q{z2A -x) + exp{-KAZ2A)Q{x - Z2A)]dx 

Gba^Fb r Q{z^b - x){z - xf 
Jo 

[exp{-KBx)Q{z2B - x)-\- exp{-KBZ2B)^{x - Z2B)]dx 

These integrals can be easily taken in an analytical form. 
For Gaa one can get 

• For z > ZsA with arbitrary Z2a 

3 

Gaa = XI ^^P^i 

i=0 



where 



(-1)*3! /-^sA . 

PAi = -T77Z TTtFa / X'[exp{-KAZ2A)Q{X - Z2a) + 

i\\6 — ^)! Jo 



;-i)^3! 

i!(3-i)! 

ex^{—KAx)Q{z2A — x)]dx 



• For z < min{z2A: Z3a} 

Gaa — Fa I {z — x)^ eyo^i—KAX^dx 

Jo 

has been already calculated. 
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For Z2A< Z < Z-iA 



Gaa = J2 ^'PAi + Fa— — exp{-KAZ2A) 



where 



PAi = FA—- -7 / x'exp{-KAx)dx 

il(6 — in Jo 



!(3 - i) 

All pAi are constants and can be calculated analytically. 
For Gb4: one can get 

• For z > ZsB with arbitrary Z2b 

3 
i=0 



where 



(-1)'3' f^SB . 

PBi = TTT^ TTtFb / X'[ex.p{-KBZ2B)Q{X - ^2b) + 

il{6 — ty. Jo 



(-1)'3! ^ /--sB^.^ 
i!(3-i)! 

exp(— irBa:)©(2;2B — x)]dx 



• For z < min{z2B, z^b} 

Gba — Fb I iz — x)^ ex.p(—KBx)dx 
Jo 

has been already calculated. 

• For Z2B < Z < Z3B 

^ i , n ~ Z2b)^ , N 

Gba = }_^zpBi + Fb exp(-AB^2s) 

i=o ^ 

where 

3!(-l)' /•^2s 



3!(— i) /"^zs 

Psi = Fb^tt^ ^ / x' fi-icp{-KBx)dx 

l\{6 — l)\ Jo 
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All psi are constants and can be calculated analytically. 

One can easily see that z^a < Z2a, z^b < Z2b and the third opportunity 
can not take place. 

One can calculate by the same way all further iterations^. 

We see that due to the difference of the relations between parameters 
ZAi, zbi even the functional form of the iterations will be different. But it is 
quite easy to see that calculations of further iterations will be done absolutely 
analogously. Moreover, the functional forms of all odd further iterations will 
be similar. The functional forms of all even iterations will be also similar. 

Then one can extract these functional forms and it is necessary to obtain 
the values of parameters ZAi, zb-i in these forms. 

Moreover, one can put one and the same values of "final" parameters za, 
zb in equations for Ga, Gb, Oa, Qb- This leads to the following result 

In the initial expressions 

Ga = Fa Hz - exp(-r^(G^ + GB))eAdx 
Jo 

Gb = Fb I\z - x)^ exp(-r^(G^ + Gb))^^^^^ 
Jo 

Then ^ 

Ga = Fa Q(za — x)(z — x)^ exp(—KAx)dx 
Jo 

Gb = Fb / ^{zb — x){z — x)^ exp{~KBx)dx 
Jo 

These integrals can be easily taken in analytical form. 
For Ga one can get 

• For z > Za 

3 

Ga = Y1 ^'PAi 
i=0 

where 

(-1)*3! /"^A . 

PAi = ^777: t^Fa x' exp{-KAx)dx 

i\{3 — ly. Jo 

■^That's why we left all possibilities in the last expressions. 
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• For z < za 

Gaa = Fa I (z — x)^ exp(—KAx)dx 
Jo 

has been already calculated. 

All pAi are constants and can be calculated analytically. 
For Gb one can get 

• For z > zb 

3 

Gb = J2 ^'Pbi 

where 

(-1)*3! f^B . 
PBi = -T77y TT^-B / a;*exp(-K'Ba;)(ia; 

• For z < Zb 

Gb4 = Fb (z — x)'^ exp(—KBx)dx 
Jo 

has been already calculated. 

All PBi are constants and can be calculated analytically. 

As the result now we know the functional form for Ga, Gb, Oa, Sb- The 
condensations equations system now determines the characteristic values za, 
Zb- 

The step-like behavior of occurs due to the avalanche character of the 
vapor consumption. The avalanche character of the vapor consumption has to 
be considered as the main feature of the first order phase transition kinetics. 
This feature was described in in details. Meanwhile it is possible to make 
the theory more accurate taking into account the regime of the droplets 
growth. 

Note that all functions if are lying between exp(— x^) and exp(— x'^) after 
the suitable renormalization. Namely, they are more sharp than exp(— x'^) 
and more smooth than exp(— x^). Note that 

/•oo rco 

/ exp(— x^)(ix ~ / exp(— x^)(ix ~ 0.9 = D 
Jo Jo 

Then one can make results of the theory more accurate if we take 

• instead 0(zjA ~ the function Q{DziA — x) 
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• instead Q{ziB — x) the function Q{DziB — x) 

• instead Q{x — zia) the function DQ{x — Dz-ia) 

• instead 9(x — Zis) the function DQ{x — Dzib) 

and multiply all parameters ZiA, ZiB on D. 

We can see that the generalization of the presented approach to the mul- 
ticomponent case is quite evident. But here the number of unknown parame- 
ters za will be equal to the number of the sorts of heterogeneous centers. The 
problem how to solve this system is rather actual and here we shall propose 
a way to do it. 

We shall choose the sorts of heterogeneous centers as to have 

Ta>Tb>Tc>Td> .... 

Then 

Za< Zb < Zc < ZD < 

At first we shall determine za- Note that the behavior of the system after 
Za- Then we can put all all F^, Tc etc. to be equal to Va and zb, zc etc. to 
be equal za- Then we have only one parameter za and can determine it from 
one algebraic equation analogous to those described in ||^ (the situation is 
absolutely analogous to the intermediate situation). 

Now Za is determined and Ga, Oa are known. It allows to go to the 
determination of zb- It can be done by the same procedure. We have to 
consider Ga, Oa as known values and put all zc, zd etc. to be equal to zb- 
Then the analogous equation can be solved and it gives zb, Gb, Ob- This 
procedure can be continued. Note that at first this way was described in p 
for the special monodisperse approximation. 

3 Special monodisperse approximation 

The task to give the theoretical description of the complex systems required 
new approximations for some characteristic functions. As the result the 
monodisperse approximation for Ga and Gb was suggested in [^. The 
monodisperse approximation can be used in two variants - the first one is 
the monodisperse approximation with fixed number of droplets which is suit- 
able when Ga, Gb are really important, the second one is the monodisperse 
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approximation with a floating number of droplets involved in this approxi- 
mation. A special recipe allows to get approximation which is suitable during 
all times. 

The second variant of approximation is necessary for the systems with the 
strong hierarchy between the probabilities of droplets formation on different 
sorts of centers. Nevertheless in solution of the system by the procedure 
analogous to the already described in the end of the last section it is possible 
to use only the first type of approximation. This states the significance of 
the first type of approximation. But as far as we follow here the iteration 
procedure it will be necessary to use the second variant of approximation. 

The problems to calculate iterations appear in the second approximation: 



The justification of the monodisperse approximation has been already 
presented many times (see @], 0) and we shall use it here without discussion. 



The real advantage of consideration of the monodisperse approximation on 
the level of iterations is the possibility to calculate the error of approximation 
explicitly. This error is small. 

One can make the iteration approximation more accurate by the substi- 
tution of exp{—KAx) by 





Then 




exp(-rA(-FA + Fb)—) exp{-KAx)dx 
exp{-TB{FA + Fb)—) exp{-KBx)dx 





Jo 



and by the substitution of exp(— K^x) by 



92B = exp{-KB / exp{-TB{FA + FB)—)dx) 
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This leads to 

G2A{z)=FA{l-e2AC^))zyKA 
G2b{z) = Fb{1 - e2B{\))zVKB 

When monodisperse approximation is used it convenient to use the fol- 
lowing way to construct iterations 

GAi+i{z) = Fa I\z - xfQx.^{TA{Ux) - ^))e Ai+i{x)dx 
Jo 

GBi+i{z) ^ Fb [ {z- xfex.p{TB{Ci{x) - ^))eBi+i{x)dx 
Jo 

All other formulas remain the same. Certainly the chains of inequalities 
will be slightly violated, but this can not leads to divergence of iterations. 
Moreover, one can simply add the difference between sequential 6 to the error 
between sequential iterations. 

Now we see that in monodisperse approximation 

GiA{z)^^{l-9Ai{z/A)) 

G',b(^) = %^(1-M^/4)) 

So, we need to calculate only Oai, Obi- 

In this calculation one can use the avalanche character of the vapor con- 
sumption again. Then 

d2A{z) = exp{-KAZiA)0{z - zia) + exp{-KAz)Q{ziA - z) 

02b{z) = e-x?^{-KBZiB)0{z - zib) + exp(-XB2;)©(2;iB - z) 

where Zxai Zib are given by the previous expressions (but further ZiA, ZiB 
will have slightly another numerical values). 

Note that in practice we don't use such expression but can use the simple 
asymptotes ("essential asymptotes") 

92A = e^{-KAz) 

92B = exp{-KBz) 
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for all z. 

Really, we need these expressions only iox z < za^ z < zb- As far as the 
argument is zjA, it means that these asymptotes can not be used only when 

> or zb > 4:Za- It can take place only when Fa and Tb differ more 
than in 4^ = 64 times. The last seems to be practically unrealizable. 

The third iteration for 6 gives 

03A = exp{-KA / exp(-rA(FA(l - exp{-KA-) 
Jo 4 

©(-| + Zia) - exp{-KAZiA)0{-ZiA + |)) + 
FB{l-eM-KB^M~ + z,B)- 

X 

exp{-KBZiB)Q{-ZiB + j))))dx) 

OzB = exp(-XB / exp(-rB(FA(l - ex^{-KA-) 
Jo 4 

+ Zia) - bM-^aZiaM-Zia + |)) + 
FB{l-exp{-KB^)e{~ + z^B)- 
exp{-KBZiB)0{-ZiB + ■^))))dx) 



One can directly see that 

fz/4 
10 



Ka 



GAi+i{z) = Faz^ r eM-'^AiGAiix) + GBiix))) 
Jo 

ex.p{-KA [ ex.p{-TA{GiA{x') + GiB{x')))dx')dx = 
Jo 

[I - exp(-ii'^ / exp(-r^(G'j^(x) + GiB{x)))dx)) 
Jo 



GBi+i{z) = Fbz^ r^\M-^B(GAi(x)+GBi(x))) 
Jo 

rx 

exp{-KB / exp{-TB{GiA{x') + GiB{x')))dx')dx = 
Jo 

^ (1 - exp(-XB / exp{-VB{GiA{x) + GiB{x)))dx)) 
Jo 



K 
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Then Gas, Gbs have the same structure as Ga2, Gb2 have. Certainly the 
numerical values of parameters zai, zbi can be slightly another. All other 
iteration approximations will have the same form. 

Now we can reformulate the system of the condensation equations because 
we know the functional form of Ga, Gb- It will be the following 

Ga{z) = ^(1 - eM-KA / exp(-r^(G'^(a;) + GB{x)))dx)) 
Ka Jo 



Ka 

Gb{z 



^(1 - eM-KB r'\M-^B{GA{x) + GB{x)))dx)) 
Kb Jo 



In the avalanche approximation it will be the following 

FaZ^ 

Ga{z) = ^(1 - exp(-i^^^/4)0(^^ - z/A) - exp(-i^^^^e(^/4 - za)) 
Ka 

Gb{z) = ^(1 - cxy>{-Kbz/A)Q{zb - z/A) - exv{-KBZBQ{z/ A - zb)) 
Kb 

where za, Zb are defined by the previous relations. 
In the approximation of essential asymptotes 

FaZ^ 

GA{z)^^{l-eM-KAz/A)) 
Ka 

FnZ^ 

GB{z) = ^{l-eM-KBzlA)) 
Kb 

These algebraic equations can be easily solved. 

Now one can mention some ways to go away from the avalanche approxi- 
mation. In the approximation of essential asymptotes everything is clear. In 
the avalanche approximation one can consider the number of the free hetero- 
geneous centers as the smooth function of z and approximately consider 

Ga{z) = ^(1 - eM-KA r^\xp{TA{GA{x) + GB{x)))dx)) 
Ka Jo 

Gb{z) = ^(1 - eM-KB r^\xp{TB{GA{x) + GB{x)))dx)) 
Ka Jo 

Then Ga ~ const z^, Gb ^ constz^ , we have the monodisperse approxi- 
mation with fixed number of droplets and need only to calculate the integrals 
of the type j'^ ex]){—x^)dx. This leads to the necessity to take 
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• instead Q{ziA — x) the function Q{DziA — x) 

• instead Q{ziB — x) the function Q{DziB — x) 

• instead 6(x — ZiA) the function DQ{x — Dzi^) 

• instead Q{x — zib) the function DQ{x — Dzib) 
and to multiply all parameters Zia-, ZiB on D. Here 



The generalization of the presented approach to the multicomponent case 
is quite evident. Our actions are absolutely identical to those described for 
the avalanche model. The only thing to stress is that at every step we can 
use "essential asymptote" for all sorts of heterogeneous centers whose Ga 
are still undetermined. Those Ga which are already determined have to be 
taken explicitly. Note that at first this way was described in 0. 

As the result we have to conclude that approximations of the avalanche 
consumption or the special monodisperse approximation lead to the simi- 
larity of the iterations of the high order (at least in their functional form). 
So, we see now the correspondence between the iteration approach and ap- 
proximations used here. The numerical calculations show that the errors of 
these approximations are negligible in frames of the accuracy of the modern 
experiment. 



References 

[1] Kuni F.M., Grinin A.P., KoUoidn. journ. vol.46 (1984), p.460 

[2] Wakeshima H., Time lag in self nucleation J.Chem.Phys., 1954, v. 22. 
N.9, p.1614-1615 

[3] F.M. Kuni, A.P.Grinin, A.S. Kabanov, KoUoidn. journ. (USSR) V.46 
(1984), p.440 

[4] V.B.Kurasov, Deponed in VINITI 2594b95 from 19.09.1995, 28p. 




20 



[5] Kurasov V.B., Universality in kinetics of the first order phase transitions, 
SPb, 1997, 400p. 

[6] Kurasov V.B., Development of the universality conception in the first 
order phase transitions, SPb, 1998, 125p. 



21 



